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ABSTRACT 



There exist both continuum and lattice regularizations of gauge theo- 
ries with fermions which preserve chiral U(l) invariance ("fermion num- 
ber"). Such regularizations necessarily break gauge invariance but, in a 
covariant gauge, one recovers gauge invariance to all orders in perturba- 
tion theory by including suitable counter terms. 

At the non-perturbative level, an apparent conflict then arises between 
the chiral U(l) symmetry of the regularized theory and the existence of 
't Hooft vertices in the renormalized theory. The only possible resolution 
of the paradox is that the chiral U(l) symmetry is broken spontaneously 
in the enlarged Hilbert space of the covariantly gauge-fixed theory. The 
corresponding Goldstone pole is unphysical. The theory must therefore be 
defined by introducing a small fermion-mass term that breaks explicitly 
the chiral U(l) invariance, and is sent to zero after the infinite- volume limit 
has been taken. Using this careful definition (and a lattice regularization) 
for the calculation of correlation functions in the one-instanton sector, we 
show that the 't Hooft vertices are recovered as expected. 



1. Introduction and Conclusion 

Every gauge theory coupled to massless fermions has an anomalous chiral current. 
Representing all fermions by left-handed Weyl fields, the Noether current associated 
with a common global U(l) rotation is classically conserved. At the one-loop level, a 
gauge-invariant definition of the current yields the Adler-Bell-Jackiw anomaly [1] 

« = ^trFF. (1.1) 

For our notation see Appendix A. The group-theoretical constant c is additive. (Each 
Weyl fermion in the fundamental representation contributes c—\. We will assume 
the gauge symmetry to be non-anomalous throughout this paper.) One can also define 
a conserved but gauge non-invariant current 

j£ = Jt-g 2 K lt , (1.2) 

where 

If a gauge-invariant regularization is used, the gauge-invariant, non-conserved current 
is defined (up to a Z factor) by the fermion bilinear 

E^i^i- (i-4) 

i 

Here % runs over all the left-handed fields. In QCD-like theories, this applies in 
particular to dimensional regularization as well as to the standard lattice regulariza- 
tion [2, 3]. 

What happens if the regulator is chiral-U(l) invariant? The U(l) current will now 
be conserved at the one-loop level. Therefore it must, when the cutoff is removed, 
coincide with the gauge non-invariant current defined in eq. (1.2). (This is true 
up to a term d u H^ u with an anti-symmetric tensor.) Since, classically, the U(l) 
Noether current is gauge invariant, this can only happen because the regularization 
itself is not gauge invariant: a chiral-U(l)-invariant regularization is, necessarily, not 
gauge invariant. 

Does this observation imply that all chiral-U(l)-invariant regularizations must 
be dismissed? To begin with, in perturbation theory the answer is no, provided 
the action contains covariant gauge-fixing (and ghost) terms. A covariant gauge 
will be assumed throughout this paper. In the presence of a longitudinal kinetic 
term, (d^A^) 2 , the theory is renormalizable by power-counting without relying on 
gauge invariance. The renormalization program reduces to an algebraic problem and 
(provided the gauge symmetry is non-anomalous) one can restore gauge invariance to 
all orders in perturbation theory by suitable counter terms (see e.g. ref. [4]). 

Beyond perturbation theory, there is an apparent conflict between chiral U(l) 
invariance of the regularized theory, and the fact that instanton-mediated amplitudes 
violate the conservation of the chiral U(l) charge [5]. It has been pointed out long 
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ago [6, 7, 8] that, in a covariant gauge, the breaking of chiral U(l) invariance can be 
spontaneous in a technical sense. The reason is that the enlarged Hilbert space of the 
gauge-fixed theory can accommodate a new Goldstone pole. The latter is unphysical 
since it originates from the part of the current j£. If the regulator is chiral-U(l) 
invariant there is, in fact, no other possibility. In this paper, we re-examine this ques- 
tion in the context of a specific lattice-regularization method. Our analysis reveals 
that a careful definition of the thermodynamical limit is necessary, just as in the 
case of conventional spontaneous symmetry breaking. This generalizes to continuum 
regularizations with chiral U(l) invariance such as, for example, momentum- cut off 
schemes (see e.g. ref. [9]) or the dimensional-reduction scheme of ref. [10]. 

The main motivation for our lattice-regularization method is that it may ulti- 
mately provide a non-perturbative definition of (anomaly-free) chiral gauge theories. 
These theories are notoriously difficult to regularize in a gauge-invariant way. In par- 
ticular, dimensional regularization is not a gauge-invariant regulator in this case. A 
gauge-invariant perturbative regularization for chiral gauge theories has been found 
only recently in the context of lattice gauge theory [11, 12]. Beyond perturbation 
theory, it is not known if non-abelian chiral gauge theories can be regularized in a 
gauge-invariant manner. (For a review of recent work in this direction see ref. [13].) 

According to the gauge-fixing approach, chiral gauge theories are defined as the 
continuum limit of a lattice theory whose action contains a covariant gauge- fixing term 
and counter terms [14, 15, 16, 17] (see also ref. [18] for a pedagogical presentation, 
and ref. [19] for a recent review on lattice chiral gauge theories). Carefully chosen 
irrelevant terms in the lattice action are essential for the existence and the continuity 
of the phase transition where the continuum limit is taken. In the abelian case, we 
showed that the lattice fermions are indeed chirally coupled to the gauge field, and 
that perturbation theory provides a valid description of the critical point [16]. 

A fully non-perturbative generalization of the gauge-fixing approach to non- 
abelian theories is not a simple task, since one has to confront the issue of Gribov 
copies. This important problem will not be addressed here (see ref. [17] for recent 
progress). Still, the method generates a systematic expansion around the classical 
vacuum and, provided the fermion spectrum is gauge-anomaly free, it provides a con- 
sistent regularization in perturbation theory. By invoking the familiar machinery of 
collective coordinates [5], it can be used to generate a systematic expansion around 
other classical solutions. In particular, it is possible to carry out an analytic calcula- 
tion in an instanton background. This allows us to address the question, first raised 
in ref. [20] , of how fermion-number violating processes are realized in the gauge-fixing 
approach. 

Let us explain the issue in more detail. The simplest chiral fermion action used 
in the gauge-fixing approach, the so-called chiral Wilson action, utilizes a set of right- 
handed spectator fields \r) one P er each left-handed field ip l L (see Appendix C for 
the precise definition). The role of the spectators is to avoid fermion doubling on the 
lattice [2, 3, 21]. Thanks to a fermion-shift symmetry they can be proven to decouple 
in the continuum limit [22]. 

Each term in the lattice action has one ip L or Xr, an d one ipL or xr field. The 
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chiral Wilson action is therefore invariant under a common U(l) rotation of a// fermion 
fields. The Noether current associated with this symmetry, 

„ Yfp L <t„ Vl + t R Xr + ■ ■ •) , (1-5) 

i 

is exactly conserved on the lattice (the dots stand for lattice terms with no continuum 
counterpart, see Appendix C). But since the lattice action is not gauge (nor BRST) 
invariant in the gauge-fixing approach, the conservation of J^' latt is consistent with our 
earlier general comments. Following ref. [23], we have verified through an explicit one- 
loop lattice calculation [24] that the current jL> latt indeed reduces in the continuum 
limit to the current defined in eq. (1.2). 

The perturbative results of refs. [23, 24] are, however, not enough to resolve 
the following puzzle, that we will refer to as the Banks paradox [20]. In short, the 
paradox stems from the fact that exact U(l) invariance implies exact conservation of 
the corresponding U(l) charge, namely, of fermion number. To be precise, consider 
a finite-volume lattice and assume that the boundary conditions respect the U(l) 
symmetry. The Ward identity that corresponds to a global U(l) rotation reads 

(SO) = q(O) = 0. (1.6) 

Here q is the fermion-number, or the U(l) charge, of O. We stress that, on a finite 
lattice, this is a rigorous result. (The invariance of the lattice measure follows trivially 
from the fact that there is an equal number of dipi and dip L Grassmann integrals, as 
well as an equal number of dxR and dx R ones.) 

The Ward identity (1.6) states that a fermion correlation function (0(xi,x 2 , ■ ■ .)) 
can be non-zero only if q — 0, namely if the number of and xr fields is equal to 
the number of ip L and Xr fields. This is true for any finite lattice spacing a and any 
finite volume, and therefore also after the continuum limit a — > and the infinite- 
volume limit have been taken. Moreover, in the continuum limit the numbers of xr 
and Xr fields must by themselves be equal, since the spectator field decouples [22]. 
Hence the numbers of ipi and ip L fields must be equal too. We have thus reached the 
paradoxical conclusion that, even though the lattice-fermion spectrum is chiral [16], 
all fermion-number violating amplitudes vanish! In other words, 't Hooft vertices [5] 
do not seem to occur. If the gauge-fixing approach would be utilized to define a 
vector-like theory such as massless QCD, the same reasoning would seem to lead 
to the erroneous conclusion that the U(l) axial charge is conserved in all physical 
processes. 

One can avoid the Banks paradox by adding to the lattice action a mass term, 
rriijip^ipL + h.c, for the physical fermions. This allows unequal numbers of ipi and 
i\) L fields to be compensated by insertions of the mass term. The broken, global U(l) 
Ward identity in a finite (lattice) volume is now 

q{0)^2m ij l^^i{x)^i{x)-h.c^ . (1.7) 
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The original lattice-fermion action corresponds to the limit m — > 0, where m denotes 
generically the magnitude of the mass terms. We see that fermion-number or axial- 
charge violating amplitudes can be non-zero provided they behave like m/m in the 
limit m — > 0. Now, we would not expect an "m/m" behavior in a finite volume. 
(Being valid for m = 0, eq. (1.6) in fact implies that this is impossible in the presence 
of a lattice cutoff.) The conclusion is that, in order to reproduce correctly 't Hooft 
vertices, the infinite-volume limit and the m — > limit must not commute. Hence, 
the U(l) lattice symmetry must be broken spontaneously. 

Anticipating spontaneous symmetry breaking (SSB) of the U(l) lattice symmetry 
we define the thermodynamical limit as the infinite- volume limit, followed by the limit 
m — > 0. In this limit the momentum-space U(l) Ward identity reads, for any p^ ^ 0, 

i P »(ji Mt (p)0)=q(0) . (1.8) 

Here J^ ,latt (p) is the Fourier transform of J^' latt (a;). Unlike eq. (1.7), because here 
Pn 7^ 0, the explicit m-dependent term now vanishes for m — > 0. This is explained 
in more detail in the discussion section. The identity (1.8) holds in the regularized 
theory as well as in the continuum limit. (On the lattice, the p^ factor on the left- 
hand side is modified by 0(ap 2 ) terms. Again, remember that there is no room for 
an anomalous term since the current J^ latt is exactly conserved on the lattice.) 

Since 't Hooft vertices do exist, this means that there are operators for which 
the right-hand side of eq. (1.8) is non-zero. Hence the left-hand side must contain 
a Goldstone pole. As explained earlier, in the continuum limit the Goldstone pole 
comes from the part of the conserved U(l) current, and is unphysical [6, 7, 8]. 

In this paper we calculate instanton-sector fermion correlation functions on the 
lattice, in the semi-classical approximation. We start from a lattice-fermion action 
with an additional, small, mass term that breaks explicitly the unphysical U(l) sym- 
metry. Taking the infinite- volume and continuum limits, followed by the m — > limit, 
we show that the anticipated 't Hooft vertices are recovered. The "m/m" nature of 
fermion-number and axial-charge violating amplitudes is manifest in our calculation. 

The paper is organized as follows. In order to minimize technicalities we begin 
in Sect. 2 with one-flavor QCD where, instead of the usual gauge- invariant lattice 
definition, we define the theory via the gauge-fixing approach in the special case that 
the (left-handed) fermion spectrum happen to contain one field in the fundamen- 
tal representation and one field in the complex conjugate one. In Sect. 3 we work 
out the anomaly-free SO(10) theory as the prototype of a truly chiral gauge theory. 
Our conclusions are summarized and discussed in Sect. 4. In particular, we show 
that the phase of the 't Hooft vertex follows the phase of the applied mass term, 
thus demonstrating explicitly the existence of the continuously degenerate ground 
states associated with SSB. Notations are listed in Appendix A, elements of SO(10) 
group theory are discussed in Appendix B, and lattice definitions are collected in 
Appendix C. The construction of propagators in the presence of approximate zero 
modes is discussed in Appendix D. 
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2. One-flavor QCD using the gauge-fixing approach 

We begin with the simple example of one-flavor massless QCD, an SU(N) gauge 
theory coupled to one Dirac fermion in the fundamental representation. The anomal- 
ous current of eq. (1.1) is in this case the axial current. Let us recall what the 't Hooft 
interaction of this theory is. In a fixed instanton background, the massless contin- 
uum Dirac operator p has one left-handed zero mode u(x) = Plu(x). Therefore the 
Weyl fields and ip R each have a zero mode. (In an anti-instanton background, 
the Weyl fields with zero modes are ipR and ip L .) The basic axial-symmetry violating 
correlation function is 

(Mx)^n(y)) = u(x)u\y) Bet', (2.1) 

where the expectation value denotes Grassmann integration only, and Det' is the 
(renormalized) fermion determinant with the zero mode removed. Our objective will 
be to recover this result starting from a lattice action with exact axial U(l) invariance. 

The remaining step in a complete semi-classical calculation is the integration over 
the gauge and ghost fields. This raises no new conceptual issues and therefore we will 
skip the details. We recall that the integration (or lattice-sum) over the instanton 
position recovers momentum conservation. The gaussian integration over the non- 
zero gauge, fermion and ghost fluctuations leads to the replacement of the lattice's 
bare instanton action, 8n 2 /g 2 , by the renormalized one, 8n 2 /g 2 (p), where p is the 
instanton's size. 

One-flavor QCD has a gauge-invariant lattice definition. When using ordinary 
Wilson fermions, the lattice fermion action is not invariant under axial transforma- 
tions, and the paradox described in the introduction does not arise. (In the continuum 
limit one reproduces the axial anomaly [3], while non-singlet axial symmetries are re- 
covered [25].) In principle, it should be possible to define one-flavor QCD using the 
gauge-fixing approach, too. While this has many disadvantages compared to the 
gauge-invariant definition, it has the interesting property that the paradox described 
in the introduction occurs. By working out one-flavor QCD we are able to address, 
with minimal technicalities, the main issue of this paper — how a global symmetry 
of the lattice path integral can be broken in the continuum limit. 

The lattice construction of the chiral Wilson action begins with enumerating the 
left-handed fields of the target theory. For one-flavor QCD we have two Weyl fields 
ipL and ip c L in the fundamental and anti-fundamental representations respectively. As 
explained in the introduction, one also needs two right-handed spectator fields that 
decouple in the continuum limit. These may be denoted xr and Xr- m the case at 
hand, we may take advantage of the Dirac nature of the target theory, and trade the 
left-handed anti-fundamental field with a right-handed fundamental one ipl — > ip R , 
ip c L — > i/jr. With a similar tradeoff for the corresponding spectator field, the lattice 
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fermion action density can be written in the following matrix form 



Xr i>R Xl 
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\ Xr J 



(2.2) 



We use hats to denote lattice derivatives. (For the precise definitions of lattice deriva- 
tives and currents see Appendix C.) Observe that the middle two-by-two block in the 
above matrix operator resembles the massive continuum Dirac operator. 

For orientation, we recall that in the conventional definition of Wilson fermions 
there are of course no spectator fields, and covariant Wilson terms are placed in the 
same block-entries as the mass terms in the above expression. This removes the 
doublers in a gauge invariant way, while axial symmetry is lost. 

In order to later accommodate truly chiral gauge theories, the doublers are re- 
moved here by introducing spectator fields and coupling them to the original fermions 
via the (free) lattice laplacian □. Since the Wilson terms now couple fields with differ- 
ent gauge-transformation properties, they lead to breakdown of gauge invariance [26]. 
In lattice perturbation theory gauge invariance is regained by adding suitable counter 
terms, and the renormalized diagrams describe one interacting Dirac field, the quark, 
and one free Dirac field, the spectator. As is usually the case for symmetries broken 
by the lattice regularization, the above is true provided the external momenta are 
vanishingly small in lattice units. The choice of a free lattice laplacian in eq. (2.2) 
implies the shift symmetry of the spectator field [22], which reduces considerably the 
number of counter terms. In particular, there are no counter terms of the form ip R xl, 
etc., and the spectator field decouples in the continuum limit. 

Let us now examine the U(l) symmetries of the action (2.2). Dropping both 
the Wilson and the mass terms, the action would be invariant under four separate 
U(l)'s — a fermion- number symmetry for each Weyl field. With the Wilson terms 
in place, the action is still invariant under two U(l)'s. Finally, for m ^ 0, only the 
invariance under a common U(l) rotation is left. This invariance corresponds to the 
baryon-number symmetry of QCD. 

The additional U(l) symmetry at m = transforms Xl and ipn with (say) charge 
+1, and xr an d with charge —1. This is the chiral symmetry that leads to the 
Banks paradox [20]. The existence of this lattice symmetry would seem to lead to 
the (erroneous) conclusion that the axial charge is conserved in massless (one-flavor) 
QCD. As mentioned in the introduction, refs. [23, 24] already showed that the anomaly 
appears in the triangle diagram as expected. But this still does not explain how the 
axial charge is not conserved in physical processes. 

We will now answer this question through an explicit calculation. We calculate 
the lattice-fermion two point function in the semi-classical approximation for (small) 
m > 0. Taking the infinite-volume and continuum limits, and finally sending m — > 0, 
we find that the 't Hooft interaction (2.1) is recovered. 

As discussed in the literature [6, 7, 8], the continuous degeneracy of ground states 
associated with SSB of the axial U(l) is parametrized by the vacuum #-angle. The 



chiral Wilson action (C.I) is invariant under a CP transformation, and this remains 
true in the presence of the mass term introduced in eq. (2.2) above. Therefore the 
calculation in this section (as well as in section 3) corresponds to a vacuum angle 
9 — 0. The case of general 6*-angle is explained in the discussion section. 

We start with a continuum regular-gauge instanton field A^(x) whose size p is 
very large in lattice units, p ^> a. (Singular-gauge instantons are suppressed in the 
gauge-fixing approach by the irrelevant terms in the lattice action, see Appendix C.) 
The lattice gauge field may be defined as U^x) = exp(iaA^(x)). This is a smooth 
configuration. (By this we mean that U^(x) — I = 0(a/p) and U^(x) — U^(x + u) = 
0((a/ p) 2 ).) For this lattice gauge field and fixed m > 0, we denote the matrix 
operator in eq. (2.2) by 75.M. Note that according to this definition Ai is hermitian. 

In the formal continuum limit, the lattice operator in eq. (2.2) goes over to a 
continuum Dirac operator D(m) = ^y 5 H(m) which depends on the original instanton 
field Afj,(x). One obtains D(m) by dropping the (irrelevant) Wilson terms and re- 
placing the lattice difference operators and by the corresponding continuum 
derivatives. For m > 0, D(m) describes a massive quark (made of iPr,l) whose Dirac 
operator is p+m, and a decoupled, free massless spectator field (made of Xr,l)- Thus, 
D{m) has no zero modes. We will denote by G(m) the propagator of the hermitian 
operator H(m). Both 1 and G(m) admit a standard spectral decomposition. 

Let us now be more precise about how the physical matrix element is obtained. 
One has to multiply the correlation function in eq. (2.1) by (normalized) wave- 
functions fl(x) and f2(y) and integrate (or sum) over x and y. The physical observable 
is the gauge-field functional average of Det(.M) A^ _1 7s 1/2)- We will denote the 
generic virtuality of the external legs by Q 2 . As explained below, the matrix element 
is dominated by instantons of size p 2 ~ Q~ 2 ^> a 2 , where the last inequality follows 
because in the continuum limit a 2 Q 2 — > 0. 

Since the fermion determinant will be 0(m), we are interested only in the 0(l/m) 
piece of the propagator (s). We claim that 

lim (A| M-H |/ 2 )^ ar = (A| G(m) 75 |/ 2 )Ztu m ■ (2-3) 

Here "singular" denotes the 0{l /m) piece. The reason why eq. (2.3) is true is that 
a non-zero m affects only eigenfunctions with eigenvalues A in the region |A|<m. 
Indeed, for all (lattice or continuum) eigenfunctions with, say, A 2 > Q 2 , the effect of 
m > will be bounded by m 2 /Q 2 to some positive power. Therefore they do not 
contribute to the 0(1/ m) term. For A 2 < Q 2 , the difference between each continuum 
eigenfunction and the corresponding lattice eigenfunction is bounded by a 2 Q 2 to some 
positive power. Since m > 0, the inverse eigenvalues A" 1 are bounded from above, 
and the contribution of the entire low-energy lattice spectrum approaches smoothly 
the continuum one for a — > 0. (For the pairing of the lattice and the continuum 
eigenfunctions we may momentarily assume a very large, but finite, volume, thus 
making the spectrum discrete; alternatively eq. (2.3) can also be justified directly in 
the infinite- volume limit.) 
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In the infinite- volume limit, the massive continuum propagator satisfies 



G(x,y;m)-f 5 



-u{x)u\y) + 0{l). 



(2.4) 



Here we show explicitly only the term that diverges for m — > (compare eq. (D.15)). 
With the understanding that the matrix element is to be taken between smooth 
wave-functions f± t 2 as described above, we thus have 



For the fermionic determinant similar arguments lead, after renormalization, to 



The explicit factor of m again comes from the (approximate) zero mode, while the 
0{m) terms account for the change in the continuous spectrum due to m. Putting 
this together we thus obtain 



Finally, taking the limit m — > we recover eq. (2.1). Eq. (2.7) reveals the "to/to" 
nature of the 't Hooft interaction. 

The familiar eq. (2.4) above is particularly simple, and has been invoked primarily 
for pedagogical reasons. In Appendix D we show how to handle perturbations that 
lift the zero modes, but are not proportional to the identity matrix and/or are not 
spatially constant. This more general formalism will be necessary in the next section. 
For a few more details on the calculation of the determinant see Comment 2 in 
Appendix D. 

Our instanton calculation was done in the semi-classical approximation, as is 
routine in the continuum. Since we have somewhat expanded its scope by using a 
specific lattice-regularization method, we will briefly review the justification for the 
semi-classical approximation. 

Consider a fermion- number violating amplitude with only a minimal number of 
fermions, and no other particles, as the external legs. (Here we wish to avoid the 
controversy about whether the fermion- number violating cross section could become 
large at very high energies due to multi-boson final states.) As before, denote the 
generic virtualities of the external legs by Q 2 . For instanton size p 2 3> Q~ 2 , the overlap 
of the zero modes with the wave-functions on the external legs will provide a strong 
damping factor. Hence the saddle point p sp of the integration over the instanton's 
size is p 2 p ~ Q~ 2 . If Q 2 is much larger than the confinement scale, the running 
coupling g r = g r {p sp ~ VQ~ 2 ) is small. This justifies the use of the one-instanton 
approximation. 

Ultimately, the most visible consequence of the anomaly in one-flavor QCD is 
that the lightest pseudo-scalar state (the "r/ meson") is not light compared to the 
confinement scale (see e.g. ref. [8]). The chiral-symmetry breaking effect obtained 




(2.5) 



lim Det(jM) = m(Det' + 0{m)) . 



(2.6) 




(2.7) 
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from the semi-classical instanton calculation is much smaller since it is controlled by 
the small parameter exp(— 87r 2 /^(p sp )). We resort to this deep euclidean regime, 
because only there are we able to apply analytic methods to accurately calculate the 
consequences of the anomaly. 

The above considerations have to do with the asymptotically-free nature of the 
Yang-Mills coupling, and therefore they are completely independent of the regular- 
ization method. Moreover, our explicit calculation has demonstrated that no uncon- 
trolled lattice artifacts occur. Finally, we note that the discretization of regular-gauge 
instantons does yield gauge-field configurations that fail to satisfy the lattice Yang- 
Mills equation of motion, but only by a small amount 0(a/p sp ). Instanton- sector 
Feynman rules that generate a systematic expansion in gf(p sp ) can be derived in the 
presence of an approximate classical solution, see e.g. ref. [27]. 

3. Chiral gauge theories 

The lesson of the previous section is that an 't Hooft vertex can be interpreted 
as an order parameter for the spontaneous breaking of the U(l) chiral symmetry in 
a regularization scheme where chiral (but not gauge) invariance is preserved. The 
introduction of a small mass term, which is sent to zero after the infinite-volume 
limit was taken, provides the necessary coupling to an "external magnetic field" and 
allows the expectation value of an 't Hooft vertex to be non-zero. This reasoning is 
valid both in the continuum and on the lattice, if one uses the gauge-fixing approach. 

The generalization of the previous calculation to 't Hooft vertices that violate 
the fermion-number symmetry of a chiral gauge theory is relatively straightforward. 
Starting from the lattice theory, a mass perturbation that lifts the fermionic zero 
modes will again allow us to keep the (approximate) zero modes under control while 
taking the infinite-volume and continuum limits. Performing next the limit m — > 0, 
we will recover the 't Hooft vertices as before. The only step which may not be 
obvious, is that a mass perturbation that lifts all zero modes exists in the continuum. 

In this section we demonstrate the existence of the necessary mass perturbation by 
working out the example of an SO(10) chiral gauge theory. (Attempting to construct 
the necessary mass perturbation for the most general anomaly-free chiral gauge theory 
may be tedious, and the SO(10) example is general enough to encompass the Standard 
Model as well as the most popular Grand Unification schemes.) In a one-generation 
SO (10) theory the Weyl fermions reside in the complex 16 representation. We 
introduce covariant derivatives (m,n = 1, . . . , 10) 



in the 16 and the 16 representations respectively. The SO(10) generators are defined 



D, 



D. 



d^ + iA^t: 



(3.1) 



via 



] = ^(l + rxOs^ + ^i-rn)!: 




N, 



MN • 



(3.2) 
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We use the 32 by 32 representation of the ten-dimensional gamma matrices given in 
Appendix B. The (continuum) lagrangian is 



C = ip L a^D fl ip L . (3.3) 

In an instanton background there are four left-handed zero modes, one for each 
quark or lepton. We will show that a suitable mass term lifts all four zero modes. To 
prepare for the introduction of the mass term we first rewrite the lagrangian in terms 
of Majorana-like fermions 

*=( €C ^ L ) , * = (^l V L eC), (3.4) 



satisfying 



^ = ^/ T C 4 C, (3.5) 



where e is the anti-symmetric two-dimensional tensor and C4 is the four-dimensional 
charge conjugation matrix (see Appendix A). The 16 by 16 matrix C, which is related 
to the ten-dimensional charge-conjugation matrix, is defined in eq. (B.3). It satisfies 
C* = C T = C^ 1 = C. In terms of the Majorana-like fields the lagrangian is rewritten 

as 

(3.6) 

where 

B„ = 7,WA + AA)=(^ "ft}. (3.7) 

Note that D\ ^ -D : unlike the QCD case, Dq is not anti-hermitian. One can show 
that 

DjC 4 C = C 4 CD* . (3.8) 

Appendices A and B contain a number of useful relations which have been used above. 
Equipped with the Majorana formulation we introduce a mass term 

f ** = y tyhWL+^LeCtf) • (3.9) 

The mass term breaks explicitly the fermion-number symmetry and, in the limit 
m — > 0, provides the "seed" for spontaneous symmetry breaking. (The mass term 
also breaks the chiral gauge invariance, see below.) The fermion operator becomes 

D(m)=D + m. (3.10) 

Eq. (3.8) holds for D(m) too. 

We will soon prove that D{m) has no zero modes, for m 7^ 0. But first, we give 
a simple physical explanation why this should be expected. Observe that 

c ® / = -i C10 r 10 , 
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where I is the two-by-two identity matrix. Introducing a 32-component spinor *&' 
whose first 16 components are equal to \l/ we may write 



y T c 4 cv = -- (y') T c 4 c 10 r 10 (i + r n )*'. 



(3.11) 



Because of Tio, the mass term can be thought of as coming from the vacuum expec- 
tation value of a Higgs field in the 10 representation. This vacuum expectation value 
breaks SO (10) down to SO (9). Since all spinor representations of SO (9) are real, 
the fermions can acquire Majorana masses consistently with SO (9) invariance. More- 
over, the sixteen-dimensional representation of SO (9) is irreducible, and therefore all 
sixteen fermions acquire a Majorana mass. 

We will now show in more detail that there are no exact zero modes for m/0, and 
that the fermion-number- violating 't Hooft interaction is recovered in the limit m — > 0. 
We describe the main steps here, relegating further technical details to Appendix B. 
In order to obtain information on the fermion propagator for m^Owe will need the 
general formalism of Appendix D, which applies to hermitian operators. We will thus 
consider the following hermitian operator (and corresponding propagator) 



V(m) 




Q(m) 




(3.12) 



Note that D{m) carries a four-component spinor index, and V{m) carries an eight- 
component spinor index. Let us first enumerate the zero modes for m — 0. There 
are the four original left-handed zero modes Ui that belong to the 16 and satisfy 
0. In addition, define 



-eCu* 



Using the left-handedness of Ui and eq. (3.8) one has 

j( \ „J„ ( N 



-C A C D 



Ui 



= 0. 



(3.13) 



(3.14) 



Therefore the Vi are left-handed zero modes of Dq that belong to 16. The propagator 
Qo(x,y) is orthogonal to all eight zero modes and satisfies 



£>o Go(x, y) = 5 A (x -y)- V(x, y) , 
where the zero-mode projector is 

/ 

Vi(x)vi{y) 



(3.15) 



V(x,y) = 



(3.16) 







V Ui (x)ul(y) J 

From these equations one can read off the relations satisfied by the chiral propagator: 

D G (x,y) = 8\x-y)-P L v i {x)v\{y)P L , 
G (x,y)D = 5\x-y)-P LUi {x)u\{y)P L . 
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The derivative acting to the left has a minus sign. 

We now turn to m ^ 0. As explained earlier, the fermion- number symmetry is 
broken explicitly. This is reflected by the fact that G(m) does not anti-commute with 
7 5 (compare eq. (3.21) below). An inspection of eqs. (3.10) and (3.12) reveals that to 
first order, the mass perturbation can have non-zero matrix elements only between a 
16 and a 16 zero-mode. Let 

\j = (vi\ m \uj) = m (vi\uj) . (3-17) 

We find, using eq. (D.15), 

G(m) = \u i )X£(v j \+0(l). (3.18) 
The first term on the right-hand side is 0(l/m). Furthermore, using eq. (3.13), 

Xij = m J uJeCuj , (3.19) 

which implies that \j is antisymmetric. In the zero-mode sector, the Majorana- 
fermion determinant is the analytic square root (pfaffian) 

det 1 ^ 2 (A) = \tijki \j Am . (3.20) 

o 

In Appendix B we prove that det (A) ^ for (almost) every embedding of the instan- 
ton in SO(10). 

Next consider correlation functions in the instanton sector. In terms of the orig- 
inal Weyl fields, 

G(m) = {**)= \ \ _ 1 V ' . (3.21) 

Let us first see what happens if we saturate two fermion fields ipL{x)ipL(y) by the 
0(l/m) part of the propagator. Using eqs. (3.13), (3.18), and noting the lower-right 
entry of eq. (3.21), we obtain a factor 

^(xK(y)A-. 1 . (3.22) 

By itself, this will give a vanishing result in the limit m — > because det 1 / 2 (A) is 
0{m 2 ). Next consider saturating the product of four fields ^L^ipLii/^L^z^^w). 
Summing over all possible contractions and paying attention to Fermi statistics we 
get 

u t (x) Uj (y)Mz>iM (A,^ 1 + A", 1 A," 1 + X^Xj, 1 ) . (3.23) 

The expression in parentheses is completely anti-symmetric in the four indices i,j, k, I. 
(If we would try to saturate 2n fermion fields for some n > 2 with the 0(1 /m) part of 
the propagator, the result would be identically zero due to the anti-symmetrization.) 
To evaluate the sum we write 

(V A « + A ^ V + ^ V) = c ^ > ( 3 - 24 ) 
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and contracting with another e^i we find 

c = I ^kiK^Kl = detV^A- 1 ) . (3.25) 
This cancels against eq. (3.20), and, in the limit m — > 0, one is left with 

(ip L (x)ip L (y)ip L (z)ip L (w)) = e ijkl Ui(x)uj(y)u k (z)ui(w) Det' , (3.26) 

which is the expected 't Hooft vertex. 

The mass term (3.9) breaks not only the unwanted fermion-number symmetry; it 
also breaks chiral gauge invariance. This, however, does not lead to any disaster; in 
fact, we know that in the UV- regulated theory the gauge symmetry is already bro- 
ken by the regulator. The crucial point is that, in a covariantly gauge-fixed theory, 
the ultra-violet behavior of the vector-boson propagator is 1/p 2 for all polarizations. 
Consequently, the theory remains renormalizable even if terms that break gauge in- 
variance are added. 

Moreover, the addition of a mass term does not change the nature of the coupling 
of the theory. The one-loop beta function is unaffected by the mass term, and so the 
theory is still asymptotically free. Also, provided we are careful to employ a mass- 
independent renormalization prescription, universality of the renormalized coupling 
should be preserved. The same considerations imply that the continuum limit of the 
lattice theory should exist for m/0, too. In the renormalized theory, a fermion-mass 
term is expected to induce a vector-boson mass term, and so unitarity will be violated 
by 0(m) effects. In the limit m — ► (keeping physical scales fixed) unitarity should 
be recovered. 

4. Discussion 

The classically-conserved chiral U(l) symmetry is not preserved by any gauge- 
invariant regularization, and in the quantized theory physical observables exist that 
violate this symmetry. In this paper we have considered an important aspect of 
regularization methods which are not gauge invariant but, instead, respect chiral 
U(l) invariance. Our concrete motivation to do so is the gauge-fixing approach to 
(chiral) lattice gauge theories. We showed that even with a chiral-U(l) invariant 
regulator, a careful treatment of the infra-red limit reproduces correctly the gauge 
invariant, chiral-symmetry-violating 't Hooft vertices. In essence, our conclusions are 
as follows. 

(I) Since the chiral U(l) symmetry is preserved by the regularization but is not 
respected by physical amplitudes, it must be broken spontaneously. (II) Therefore, 
in order to obtain the physical amplitudes, one should introduce a mass perturba- 
tion that breaks the chiral symmetry explicitly, and take the limit m — > after the 
infinite-volume limit. (Ill) The Hilbert space of the gauge-fixed theory will contain 
a corresponding Goldstone pole, but it is unphysical, because it originates from an 
unphysical (gauge non-invariant) conserved current. 
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The second and third statements are closely tied to the first one. When there 
is SSB the thermo dynamical limit is defined by introducing an "external magnetic 
field" (here the mass term) which is switched off after the infinite-volume limit has 
been taken. As for the existence of a massless pole, it is a consequence of locality 
and the Goldstone theorem which, in the present context, has been noted before in 



The mass perturbation allows us to avoid the Banks paradox, namely the apparent 
conflict between the symmetries of the regularized theory and of the physical ampli- 
tudes. Starting from a gauge-fixed lattice theory that has the unwanted chiral sym- 
metry, we have demonstrated through explicit examples how this mechanism works 
for the anomalous axial symmetry of QCD-like theories, and for fermion-number- 
violating processes in chiral gauge theories. 

Zero modes belong to the low end of the fermion spectrum which, in general, is 
sensitive to infra-red details such as having m^fl, finite versus infinite volume, and 
choices of boundary conditions. Keeping m > provides an infra-red regularization 
for all fermionic correlation functions, and allows us to take the infinite-volume limit 
without difficulty Once in infinite volume, we obtain the correct 't Hooft vertices 
because we evidently have the correct number of (approximate) zero modes. 

As explained earlier, because of CP invariance, the vacuum angle 9 was equal to 
zero in the previous sections. Considering the one-flavor QCD example for simplicity, 
let us examine a mass term with a general axial-U(l) phase. (The generalization to 
chiral gauge theories is straightforward; see ref. [28] for a discussion of #-vacua in the 
context of the standard Wilson action.) The lattice action is now defined by replacing 
the mass term in eq. (2.2) with a new mass term pointing in an arbitrary axial-U(l) 
direction (the parameter m is real) 



By applying an axial-U(l) transformation (ip — > e~ i6,75//2 ■?/>, x ~^ e t6j5 ^ 2 x, etc.), and 
using the invariance of the lattice theory for m = 0, we can relate the value of any 
9 7^ correlation function to its value at 9 = 0. For the correlation function of 
eq. (2.1) the result is 



The subscript 9 refers to the angle of the mass term (4.1). (Note that, in eq. (4.2), 
exp(i0) may be re-expressed as exp(—iq9/2) where q = — 2 is the axial charge of 
Mr-) 

Eq. (4.2) is a rigorous result in the lattice theory. A similar relation holds in 
the thermodynamical limit. Hence the order parameter for axial-U(l) symmetry 
breaking — the 't Hooft vertex — acquires a phase equal to the phase of the applied 
mass term. This proves that the SSB ground states are indeed parametrized by the 
(vacuum) #-angle. 

The relevance of 9- vacua to the U(l) problem was discussed in detail in the litera- 
ture. In a continuum treatment the path integral can only be expanded perturbatively 
around selected classical fields, and 6*-vacua must be incorporated "by hand." Here, 



refs. [6, 7]. 




(4.1) 




(4.2) 
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we showed that the #-angle of the 't Hooft vertex arises as an unavoidable consequence 
of the lattice regularization. 

While our explicit lattice calculation does not depend on this, let us recall some 
related observations from refs. [6, 7, 8]. Any sector with a fixed topological charge, 
such as the one-instanton sector, is a superposition of all the #-vacua. But for any 
m^O, however small, the vacuum energy density is a non-trivial function of 9. This 
explains why a unique vacuum is selected in the thermo dynamical limit. 

Had the infinite- volume limit been taken while keeping m — 0, the 9- vacua would 
have remained exactly degenerate. This prescription would yield a vanishing result 
for chiral-symmetry violating amplitudes. But since clustering would be violated, 
this prescription is inconsistent. Returning to the lattice regularization, this may be 
explained as follows. If the lattice volume is finite, V < oo, the limit m — > must be 
smooth and, hence, independent of 9. Namely, 

^r^VKoofi = ^)^0,V<oc,*=0 • ( 4 - 3 ) 

Since this must be true for all 9 simultaneously with eq. (4.2), the conclusion is that 
the 't Hooft vertex vanishes on a finite lattice, if we set m — 0. (This argument is 
really an alternative explanation of the Ward identity (1.6), but makes the role of the 
axial phase more explicit.) An implication is that, clearly, one would have to keep 
m ^ in a numerical simulation in order to recover 't Hooft vertices. 

The 't Hooft vertices are characterized by an "m/m" behavior in the thermo- 
dynamical limit. Given a chiral U(l) Ward identity, should an "m/m" behavior be 
expected from any other term except the symmetry-breaking expectation value? The 
answer is no. As a concrete example consider the following momentum-space Ward 
identity in one-flavor QCD with m ^ 0, as defined on the lattice via the gauge-fixing 
approach (Sect. 2). For \pa\ <C 1 it reads (compare eq. (C.7)) 

(J^{p) J 5 ) = "2m ( J 5 (p) J 5 ) + 2 fa) . (4.4) 

In this equation, is the conserved U(l) axial current in the limit m — > 0, and 
ipip and J 5 are the local scalar and pseudo-scalar lattice densities. As in eq. (1.8) 
the tilde denotes Fourier transform. The expectation value of ipip corresponds to the 
limit x = y in eq. (2.1). While in this limit the semi-classical calculation ceases to 
be reliable, we take the non-zero result for the 't Hooft vertex as an evidence that 

The contribution of an approximate Goldstone pole to ^JsQf) Js^J should be pro- 
portional to (p 2 + t>m) _1 where v is a dimensionful constant. The corresponding 
contribution to the Ward identity goes like m/{p 2 + vm) and vanishes in the limit 
m — > 0. The contribution of all other excitations to (^J^{p) should be less infra-red 

singular. Therefore nothing that behaves like "m/m" should arise from ^Js(p) J$}, so 
long as we are careful to keep the momentum not strictly zero. Indeed, sending p — > 
as a limit is an inherent part of the Goldstone theorem (see for example ref. [29]). 
For m->0we thus obtain 

ip»(j l ^{p)J,) = 2^) . (4.5) 
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This equation is a special case of the Ward identity (1.8). 

The problem of fermion- number violation in lattice chiral gauge theories was 
previously also addressed in ref. [30]. In the (axial) Schwinger model [31], these 
authors examined a lattice- fermion hamiltonian that has a "superfluous" U(l) global 
symmetry. They monitored the response of the fermion ground state to an adiabatic 
evolution of the (abelian) gauge field that changes the topological charge of the gauge 
vacuum. They found that a U(l) charge of the anticipated amount is produced in 
the process. 

The clash between chiral-U(l) and gauge symmetry is at the heart of their argu- 
ment. Because of the lack of exact gauge invariance at the lattice level, the initial and 
final bare vacua are not gauge-transforms of each other and their bare U(l) charges 
are different. During the evolution the bare charge is necessarily conserved. But since 
the bare charge of the ground state changes in the process, there is a corresponding 
change in the normal- ordered charge defined with respect to the ground state. 

The introduction of m ^ in our work was necessary to control the infra-red 
behavior of a dynamical gauge-fermion system that undergoes spontaneous symmetry 
breaking of a peculiar nature. In contrast, in ref. [30] only the response of the spectrum 
of the axial Dirac operator to an external gauge field was considered, and so it was 
not necessary to introduce the mass perturbation. 

In conclusion, in this paper we have demonstrated convincingly that, in spite of 
the exact chiral U(l) invariance of the lattice action in the gauge-fixing approach, 
fermion-number violating processes do occur, thus resolving the questions raised in 
ref. [20]. 
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Appendix A. Notation 

The group generators are normalized according to 



tr T a T b 




(A.l) 



The dual tensor is 




(A.2) 



The topological charge is 




(A.3) 
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The hermitian gamma matrices obey the Dirac algebra 

{7 M , lu) = V , tJ = 7^ • (A.4) 
In four dimensions we use the representation 

(j M = (<r, i) , ct m = <rj = (<r, -i) . (A.6) 

The chiral projectors are 

P fl = (l + 7 5 )/2, P L = (l- 75 )/2. (A.7) 

Charge conjugation matrices play a key role in the Majorana formulation of Sect. 3. 
In any even dimension the charge conjugation matrix is defined by (see e.g. ref. [32]) 

C% = -llC , (A.8) 

and satisfies C^ 1 = = C T . In 8n + 2 and 8n + 4 dimensions, C T = — C, while in 
8n + 6 and 8n dimensions, C T = C. For the above four-dimensional gamma-matrices 
the charge-conjugation matrix can be chosen as 

C, = 737i = ( J ° ) • (A.9) 
It is unique up to a sign. The two-dimensional anti-symmetric tensor (with e±2 = 1) 

c = ^2=^_ 1 J J • (A.10) 
For any even dimension and /i^ one has 

C%7, = -7JC-7, = 7?7^ = (l^fC = -(w) T C . (A.11) 
In four dimensions one has 

C 7m (1±7 5 ) = [^(1T75)] T , (A.12) 
a relation which generalizes to 4n dimensions. 



Appendix B. SO(10)-ology 

We define the ten-dimensional gamma matrices by the following tensor products 

I\ = a-i ® ax ® ax ® <7i ® <7i , 
T 2 = a 2 ® G\ ® G\ ® G\ ® G\ , 
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The ten-dimensional charge conjugation matrix is 

C 10 = iC®a 2 , (B.2) 

where 

C = a 2 <g 03 <E> 02 <8> 03 • (B.3) 
Notice that C10 anti-commutes with Tn. With the above definitions, eq. (A. 11) reads 




t MN J ' \ 




(B.4) 



Equivalent ly 



CS MJV — T, MN C — T, MN C, (B.5) 

a relation which is needed for the derivation of eq. (3.8). 

If D satisfies eq. (3.8) and has no zero modes, its inverse satisfies 

G(x,y) = (*(x)V(y)) = -C 4 CG T 0/,:r)C4C. (B.6) 

By taking a suitable limit, this generalizes to the case where D has exact zero modes 
and G(x, y) is constructed from the non-zero modes only. 

We now show that, for almost every embedding of an instanton in SO(10), 
det (A) 7^ (cf. eq. (3.19)), i.e. the mass term of Sect. 3 lifts the four zero modes. 
We will first show that det (A) 7^ for a particular embedding. We introduce 16 by 
16 matrices <Sj(Z), % — 1, . . . , 3, I — 1, . . . , 4, defined to be tensor products of four two- 
by-two matrices with <7j as the /-th factor and the identity for the rest. The SO(10) 
generators S MAr with m,n = 1, . . . , 4, generate two SU(2) groups. For <S 3 (2) = ±1 we 
label them SU(2) i?)i . 

Each zero mode is written explicitly as Ui = Ua^p^p^x) , where a, Pi, /3 2 , fa, At = 
1,2. Here a is the spin index. The other four indices correspond to the tensor prod- 
uct that defines the SO(10) generators in the 16 representation. Assuming that the 
instanton resides in e.g. the SU(2)^ subgroup defined above, the zero modes have 
<S 3 (2) = —1 (equivalently (5 2 = 2). Their SU(2) L index is j3\. Using the explicit 
solution of the isospin one-half zero mode (for a regular-gauge instanton) we have 

U a /3il32/33l34,i( x ) = N <Wi<W CW 4 ,i f(x 2 ) , (B.7) 
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where f(x 2 ) = (x 2 + p 2 )~ 3 / 2 and J\f is a normalization factor. The constants Op 3 p 4j i 
define the four independent zero modes. We will label them by the eigenvalues of 
5 2 (3) and <S 3 (4). Replacing the index i — 1, . . . , 4, by a pair of indices Ti, r 2 = 1, 2, 
we take Oj3 3 p^ T1T2 = p-sni.o'z) i3 4 t 2 - I* 1 matrix notation, O = ia 2 ® 03, and O 
has similar properties to the four- dimensional charge-conjugation matrix. Putting 
together eqs. (3.19), (B.3) and (B.7) we get 

\ T1T2 , T , T , = -mj\f 2 { d 4 x f 2 (x 2 ) tr (e 2 ) (O r (ia 2 0a 3 )O) , , 

= m(ia 2 ) T1T [(a 3 ) T2 ^, (B.8) 

which proves det (A) 7^ for this special case. On the first row, the explicit ia 2 <S> 03 
comes from the last two factors in the tensor product (B.3), while tre 2 comes from 
the first factor in this tensor product, and the explicit e in eq. (3.19). The transition 
from the first to the second row implicitly defines the normalization constant. 

Suppose now that a global rotation Re SO(10) is applied to the above special 
embedding of the instanton. The new zero modes are u[ = Rui. We claim that 
det (A(i?)) 7^ for almost every R. The proof is simple. Suppose on the contrary 
that det (A(i?)) = for every R in some open subset of SO(10). Since the embedding 
and, hence, det (A) are analytic functions of R, this would imply that det (A) = 
for all R. This, however, contradicts eq. (B.8) in the special case R = I. Therefore 
det (A) = may be true, at most, on a measure zero subset of SO(10). 

Appendix C. Lattice formulae 

The free symmetric lattice derivative is defined for any function f x by 

d x ,^f = (/a;+/t fx—fi) j 

where fi is a unit vector in the //-direction. The corresponding covariant derivative is 

D x ,^f — Q^^x,^ fx+fi ~ ^x-p,,n fx-fi) ? 

where U x ,^ is the link variable. The free lattice laplacian is 

□ijfl/ = — 2 (fx+fi + fx— fi 2/ x ) . 
(X 

Given a set of left-handed fields ip l L and corresponding spectator fields Xr, the chiral 
Wilson lagrangian is (compare the upper-right block in eq. (2.2)) 

£ = E (Vi ^A* + t R *A x* - y □ + ft □ xk)) > (ci) 

where r is Wilson parameter. This action is invariant under a U(l) rotation of all 
fermion fields 

^i^e-Vl, W L ^e-* a W L , Xr^^Xr- (C.2) 
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The conserved U(l) current is (see ref. [24]) 

J £i att = \ E { ^L,^u x ,^i x+fl + W L , x+fl ^uiAx 

i K 

+XR, x VnXR, x+fl + Xr, x +^hXr, x (C.3) 

~~ r {i )l L,xX l R^x+(i + XR^L^+p, ~ ^L,x+/iXfl,x ~ XR,x+fi' l l J L,x) j ■ 

It satisfies the conservation equation 

E(^r-o=°- ( c - 4 ) 

In the special case of one-flavor QCD let us introduce Dirac fermions ip = (iPr,iPl), 
i> = (tPl^r), X = (Xr, Xl), X = (Xl,Xr)- The axial transformation is 

^-e"*^, x^e^x, ?^?e-^ 5 , x - X^ h , (C.5) 

and the axial current is (note that in eq. (2.2) we set r = 1) 

J l% = \^x^U x ^ x+ii + ^ x+ ^ 5 ^Ul^ x 

-Xxl5lt^Xx+f, ~ Xx+^lfiXx (C.6) 
-r(lp x l5Xx+il ~ Xxl^x+iX ~ i'x+^Xx + Xx+fi^x 

For m/0 the axial current satisfies the partial-conservation equation 

\ E( JJS. - J™m) = - 2mJ ^ ■ ( c - 7 ) 

a n 

The difference operator on the left-hand side (the free backward derivative) becomes 
(1 — exp(— iap^))/a = ip^ + . . ., in momentum space. The local scalar and pseudo- 
scalar lattice densities are ip x ip x and J x ^ = ip^^ipx. As usual they are related by 
an axial rotation. They do not mix with the corresponding spectator-field densities 
thanks to the shift symmetry [22]. 

We now explain why singular-gauge instantons are suppressed on the lattice by 
the gauge-fixing action of ref. [33]. Recall that, in a singular gauge, the instanton's 
vector potential near the gauge singularity (located at x — 0) is ~ ^(x)d^ (x) / g 
where $(x) = a^x^/\x\. The magnitude of this vector potential grows like l/(p|x|). 

On the lattice let us make the (bare- field) rescaling A XtfJj — > A' = goA x ^. Drop- 
ping the prime, we expand the link variable as U Xjtt = exp(iaA XtfJj ) . The gauge-fixing 
action contains the expected longitudinal kinetic term 

i tr (?^) 2 ^?(? a ^) 2 - (a8 » 
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plus irrelevant terms. Here £o is the bare gauge-fixing parameter and index summa- 
tions have been shown explicitly. The leading irrelevant term that contributes to the 
classical potential is 

iH?^?4 (a9) 

The irrelevant terms break BRST invariance, and so there is no reason that regular- 
gauge and singular-gauge instantons will have the same lattice action. 

Consider now some lattice discretization of the singular-gauge vector potential. 
Inevitably, the (rescaled) vector potential will be 0(1/ a) in the hypercube(s) contain- 
ing the point x = and the vicinity. For such a vector potential U Xjfl — I = 0(1). 
The positivity of expression (C.9) and (since there are infinitely-many other irrelevant 
terms) of the gauge- fixing action as a whole [33], guarantees that the lattice action 
will be an 0(1) quantity times 1/qq. In the continuum limit go — > any such lattice 
configuration is suppressed. 



D. Continuum propagators in the presence of approximate zero modes 

Let H = Hq + aV be a Schrodinger-like (i.e. elliptic and self-adjoint) operator 
in a d-dimensional open infinite space. We assume that H has a finite number of 
(normalized) zero modes Ui(x), and that H has no zero modes. The inverse of H, 
denoted G, is defined by 

HG = 1, (D.l) 

where both sides are considered as operators acting on a suitable Hilbert space. In this 
appendix we explain how to construct a systematic approximation for G. (Eq. (D.l) 
may be rewritten in the familiar form HG(x,y) = 5(x — y) by taking the matrix 
element of eq. (D.l) between "position eigenstates" (x\ and \y).) 

Let us introduce the notation \i) for the zero modes of H , and \p) for the rest 
of its spectrum. (In instanton problems the zero modes are the only bound states, 
and \p) denotes the continuum of scattering states.) One has H \p) = E(p) \p). The 
propagator Go is defined by 

H G = l-\i) (i\ = \p) (p\ , (D.2) 

and has the spectral representation 

Go = \p) (Go) Pq (q\ , (G ) pq = E-\p) S(p - q) . (D.3) 

It is convenient to expand G = G(a) too using the eigenmodes of H as a complete 
orthonormal basis 

G = |<> P'j 1 (j\ + \p) G pq (q\ + \p) f pj (j\ + \j) (p| . (D.4) 

In this expansion, the basis vectors are fixed (a-independent), while the a-dependence 
is carried by the spectral functions fy, G pg and f P j. Below, we show that the spectral 
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functions can be expanded as power series in a where is 0(a), and G pq and f P j 
are 0(1). They will be used to construct approximations for G. 

We start by substituting eq. (D.4) into eq. (D.l). Using H \i) = we get 

l = aV\i)P~ 1 (j\+HG', (D.5) 

where G' consists of the last three terms on the right-hand side of eq. (D.4). By 
taking the matrix element between zero-mode states (i\ and \j) and using (p\i) = 
we get 

6 ij = aWV\k)0g + a<i\V\q)f qj . (D.6) 
By taking the matrix element between {p\ and \j) we get 

= a(p\V\k)p£ + H pq f qj , (D.7) 

where 

H pq = (p\H\q). (D.8) 

In order to solve for and f q j we have to invert H pq . Note that H pq is the continuous- 
index kernel of the operator H 1 - = \p) H pq (q\ which is, by definition, the projection 
of H onto the subspace orthogonal to the zero modes of H . The inverse of the 
projected operator is defined by (H^)' 1 = \p) (p\. The corresponding kernel 
satisfies (H^)^ 1 = \p) H^ 1 (q\. One has 

(p\ H \q) = (p\ H \q) + 0(a) = E(p) 5(p - q) + 0(a) . (D.9) 

Therefore (if- 1 ) -1 exists, and (compare eq. (D.3)) 

H' 1 = E-\p) 5(p -q) + 0(a) . (D.10) 

We are now ready to solve for and f q j. Multiplying eq. (D.7) by H'^ we get 

= -aH- 1 <p\V \k) fi£ , (D.ll) 

and substituting this in eq. (D.6) we get 

Pij = (i | (aV-a 2 V\q)H- p l (p\v) \j) . (D.12) 

Eq. (D.10) implies that H' 1 = 0(1), hence (3^ = 0(a) and f qj = 0(1). 

The parametric a-dependence of the last spectral function is determined to be 
G pq = 0(1). To show this, take the matrix element of eq. (D.5) between scattering 
states (p\ and \q), which gives 

5{p -q) = H pp ,G p , q + a(p\V \k) f* qk . (D.13) 

Since we already know that f q j = 0(1), the last term on the right-hand side is sub- 
leading. It follows that G pq = H~ q l + 0(a) = (Go) pq + 0(a). Moreover, by combining 
these results it follows that the spectral functions can be expanded as power series 
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in a, starting at the above-specified order for each spectral function. (See, however, 
Comment 1 below.) 

The propagator G involves fi^ 1 , and so it has a Laurent series starting at order 
l/a. To find the singular, 0(1/ a) piece of the propagator, we keep only the 0(a) 
term on the right-hand side of eq. (D.12). We get = A^i + 0(a 2 ), where 

X ijtl = a{i\V\j) . (D.14) 

This expression is recognized as the first-order energy shifts of the zero modes, as 
calculated using degenerate perturbation theory. Substituting in eq. (D.4) we obtain 

G = |i>A-.\ O'| + O(l). (D.15) 

While eq. (D.15) is all we will be using in the body of the paper, it is instructive 
to go one step further and construct also the 0(1) part of the propagator. To this 
end we approximate by (G ) pq . Eq. (D.12) then yields = A^ + 0(a 3 ) where 

X ij , 2 = (t\aV -a 2 VG V\j) , (D.16) 

which, as expected, includes also the second-order energy shifts. Next, making a 
similar approximation in eq. (D.ll) we have 

\P) U3 01 = -aGoV \i) A--\ 01 + 0(a) . (D.17) 

This involves the first-order correction to the zero mode's wave function. Putting 
everything together we find that, at 0(1), the propagator may be compactly expressed 
as 

G=|i,l)A^ 2 0',l|+Go + 0(a), (D.18) 

where 

\i,l) = (l-aG V)\i) . (D.19) 
We conclude this appendix with two technical comments. 

Comment 1. Perturbation theory is a-priori not valid, if the perturbation aV results 
in the disappearance of any bound states from the spectrum. What we have done 
amounts to showing that perturbation theory can still be used to construct an 0(a n ) 
approximation for the propagator, provided that the integrals occurring at (n + 1)- 
th order in perturbation theory converge. As already mentioned, in this paper we 
actually use only eq. (D.15). In four dimensions fermionic zero modes fall like |:r|~ 3 
or faster for |x| — > oo. Even if we perturb by a spatially constant mass term m (that 
does not vanish at infinity) eq. (D.15) gives the correct expression for the 0(l/m) 
part of the propagator. 

Comment 2. In a semi-classical calculation we also have to separate out the approxi- 
mate zero modes' contribution to the determinant. As discussed in detail in ref. [27], 
this can be done by splitting the functional integration into separate integrations 
over the amplitudes of the (approximate) zero modes, and over the orthogonal sub- 
space. To leading order, the integration over the zero modes' subspace gives rise to 
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det (Ai) (see eq. (D.14) above). The integration over the orthogonal subspace has a 
non- vanishing finite limit (after subtracting UV divergences). 
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